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1 Introduction 



One of the most precise measurements of the strong couphng constant as is provided by 
the decay rate r{Z — >■ hadrons). In the high energy hmit the quark masses may often 
be neglected. However, for precision measurements it is desirable to include also mass 
corrections of the form (m^/s)' with / = 1, 2, 3, . . .. This is particularly valid if one wants 
to predict the total cross section for e~^e~ into hadrons in an energy region where s and 
are of comparable magnitude. In the massless limit corrections up to order are 
known [jl], Terms of order a'^m'^/s, a^rn^/s 0, || and a^m^/s^ ^ are also available 
at present, providing an acceptable approximation from the high energy region down to 
intermediate energy values. For the energy region closer to the production threshold terms 
of higher order in rn^/s are necessary. In this paper a systematic approach is presented 
to compute these terms. For the moment the vector current correlator only is considered. 
The starting point is thereby the polarization function n(g^). With the method presented 
below we are able to evaluate the three-loop terms up to order (m^/g^)^ in an expansion 
of n(g2). 

The polarization function n(g^) is defined through 

{-9,.q' + q,q.) ^{q') = t J dx e^^^ {0\Tj,{x)UO)\0) (1) 

and the physical observable R{s) is related to n(g^) by 

R{s) = 127rlmn(g2 = s + ie). (2) 

It is convenient to write 

n(g2) = nW(g2) + ^^Cpn«(g2)+f^^Vn(2)(g2) + ... , (3) 

vr \ 7C J 

n(2) = ClU^2^ + CaCfU^^ + CpTmUl^^ + CfTuP , (4) 

and similarly for R{s). The colour factors {Cp = (N^ — l)/ {2Nc) and Ca = Nc) correspond 
to the Casimir operators of the fundamental and adjoint representations, respectively. 
For the numerical evaluation we set iV^ = 3. The trace normalization of the fundamental 
representation is T = 1/2. The number of light (massless) quark flavours is denoted by n^. 
In Eq. is the abelian contribution (quenched QED!) and is the non-abelian 

part specific for QCD. There are two fermionic contributions arising from double-bubble 
diagrams: For 11; the quark in the inner loop is massless, the external massive, whereas 
for n^"* both fermions have the same mass. The result for rI'^\s) is known analytically 
and will serve as check. The case where the external current couples to massless quarks 
and these via gluons to massive ones is treated in [j^ and will not be addressed here. 

The outline of the paper is as follows: In Sect, ^we will describe the procedure which 
allows the systematic expansion for n(g^) in analytic form for large external momentum 
and its implementation in a program. Sect. |^ contains the results separated according to 
the colour factors. In Sect. | the imaginary part of 11 (g^) is considered and compared to 
the result of a recent evaluation using semi-analytical methods . Finally the conclusions 
are presented in Sect. |^. 
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2 Large momentum procedure 



In this section the basic rules of the large momentum procedure are briefly described and 
their implementation in programs is considered. 

The methods which allow the expansion of Feynman diagrams containing either large 
masses (hard mass procedure) or large external momenta (large momentum procedure) 
have been used extensively in the recent past |Q . At one- and two- loop level it is in general 
still possible to perform the diagrammatic expansion "by hand" and only evaluate the 
integrals with algebraic programs. However, at three loops this is almost impossible, in 
particular if one is interested in an expansion up to high orders in ra^/s. 

Following the prescriptions for the expansion of an unrenormalized propagator 
type Feynman diagram in its large external momentum read: 

1. Generate all subdiagrams of the initial graph such that 

(a) they contain both vertices through which the large momentum enters and leaves 
the initial graph 

(b) they become one-particle-irreducible when the two vertices are identified. 

2. Taylor-expand the integrand of these subdiagrams in all small masses and external 
momenta generated by removing lines from the initial diagram. 

3. In the initial diagram, shrink the subdiagram to a point and insert the result ob- 
tained from the expansion in ^. 

4. Sum over all terms. 

The subgraphs from step ^ are denoted hard subgraphs or simply subgraphs, the reduced 
graphs, resulting from step 0, the co-subgraphs. 

For the non-planar three-loop topology all hard subgraphs are shown in Fig. |l[ The 
corresponding co-subgraphs are obtained from the full graph by shrinking the displayed 
lines to a point. The number of terms generated by the large momentum procedure 
increases rapidly with the number of loops in the initial diagram. Also, the relation 
between the expansion momenta in the subgraph and the loop momenta of the co-subgraph 
becomes non-trivial. Nevertheless, the rules listed above provide an algorithm well suited 
for the evaluation through a computer. Restricting to three-loop two-point functions with 
large external momentum and an arbitrary number of small masses, the procedure was 
implemented using the language PERL. The program can be divided into the following 
steps: 

1. (a) Generation of the relevant subgraphs and the corresponding co-subgraphs in- 

cluding the determination of their topologies, and 

(b) distribution of momenta in the subgraph and co-subgraph respecting the rela- 
tions between them. 

2. Calculation of these terms. 
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Figure 1: Possible topologies for the subgraphs of an NO-type diagram. 



In step ^ one only considers topologies, disregarding any properties of lines except 
their relative positions. Especially one neglects their momenta, masses, particle types, 
etc. For the representation of a diagram we label its vertices by integers and specify lines 
by their endpoints. A topology is thus described by the collection of its lines. 

To generate the subdiagrams we first note that the full graph is always one of the hard 
subgraphs. The remaining ones are obtained by going through the following steps: 

(i) Remove any combination of lines from the initial diagram. 

(ii) Remove the emerging isolated dots and binary vertices. 

(iii) Relabel the remaining vertices by {1, ... , ^vertices} and build all permutations. 

(iv) Compare with a table listing the basic one- and two-loop topologies, assign the 
proper topology if it matches one of the entries or, otherwise, discard the result. 

If a subgraph passes the last step, the corresponding co-subgraph is evidently fixed and 
its topology is determined in a similar way by applying steps (ii)-(iv). The result of 
this procedure is a database containing all relevant sub- and co-subgraphs, including 
information about their topology. Note that no selection criteria related to line properties 
have been applied so far. Thus, up to this point the procedure is universal and for each 
topology this part of the program has to be run only once and for all. 



In step |T^, the database resulting from [1^ is specialized to a specific diagram, including 



masses and momenta. The direction of the momentum carried by a line is connected with 
the order of the labels representing this line. Furthermore, for each subgraph the small 
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external momenta are routed in a way that they touch as few hues as possible. The output 
of step [| is then the input for step ||, together with all the necessary administrative files 
like makefiles etc. 

In step I the two FORM O] packages MATAD and MINCER are used which were 



already available and have been applied to many different problems. MINCER |]I2[ cal- 
culates massless three-loop propagator type integrals and therefore applies to the hard 
subgraphs. MATAD calculates massive tadpole diagrams up to three loops by using 



the corresponding recursion formulas obtained in and is used for the co-subgraphs 
accordingly. 

Allowing only for one internal mass scale m, the final result is then a power series 
in m^/g^, where q is the (large) external momentum. The coefficient functions contain 
numerical constants together with ln*(— g^/m^) and ln-'(— g^/yU^) with i < 3 and j < 2, 
where is the renormalization scale. 



3 Mass corrections to n(g^) 



With the method presented in the previous section the first seven terms of the expansion 
have been computed. For the three-loop case 18 initial diagrams have to be considered. 
Altogether 240 subdiagrams are produced when the large momentum procedure is applied. 
The corresponding numbers in the two- (one-) loop case are 2 (1) initial and 14 (3) 
subdiagrams. The high energy approximation for n(g^) reads in the MS scheme {Iqm = 
ln(-gVm2),/,^^ln(-gV/x2)): 
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where m is the MS renormahzed mass at scale fi^ and ( is Riemann's zeta-function with the 
values C2 = 7r^/6, Cs ~ 1.20206, (4 = vr'^/90 and C5 ~ 1.03693. B4 is a constant typical for 
massive three-loop integrals with = -f (4 - 4(2 In^ 2 + | In^ 2 + 16Li4(|) ^ -1.762800 
I3| . The overall renormalization of n(g^) is also performed in the MS scheme, i.e. in the 



expressions obtained after renormalization of m and only the poles are subtracted. The 
expressions renormahzed in the conventional QED scheme are obtained by subtracting 



n(0) given e.g. in [|T3|, || in order to obtain 11(0) 
QED can also be found in [O . 



0. The {iv? /q^Y terms in the case of 
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hadrons) 



According to Eq. (Q) the ratio R{s) is obtained by taking the imaginary part arising from 
the ln(— g^) terms of the above results. Note that starting from the quartic term logarithms 
in the mass fh appear which cannot be removed by a choice of the renormalization scale 
IX. A closer look into the method used for the calculation shows that starting from this 
order massive tadpoles appear, these being the source for such logarithms in agreement 
with the general discussion of [T^, [1^. For fi"^ = s which is the natural scale at high 



energies, the {m'^/s)^ and (m^/s)^ terms are free of logarithms. We refrain from listing 
the corresponding results which are trivially obtained from Eqs.(^-(|TO|). 

Using the relation between the MS and the on-shell mass leads to the ratio R{s) 
expressed in terms of the pole mass (Lms = ln(m^/s), L>,^ = ln(s//i^)): 
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The quartic terms are in agreement with the results of For the contribution containing 
massless fermions, -Rp"*, a comparison with the full analytical result is possible We 
found complete agreement up to the order considered. 

The results of the expansion can now be compared with those obtained via a semi- 
analytical procedure using the method of conformal mapping and Fade approximation |^ 
which are valid in the whole energy range. To conform with the conventions of [§, /U^ = Tn"^ 
has been adopted for this comparison. In Fig. ^ successively higher orders in (m^/ s) are 
included (dashed lines) and compared with the Fade result (narrow dots). By comparing 
the quadratic approximation with the one containing also terms of order {m? / sY (solid 
line) one observes that the quality of the approximation improves considerably. 

For completeness we also present in Figs, ^the results for the scale /i^ = s. This choice 
is more adequate to the high energy region and all functions i = A, NA, I, F approach 
the constants to be read off easily from Eqs.(|l3D-(16). 



Let us discuss separately the high energy and the low energy regions. For all three 
functions Ra, Rna and Ri and values between x = and x = 0.6 (x = 2m/y/s) the 
expansions including terms of order {m? / sY (or more) are in perfect agreement with 
the semi-analytical result. Conversely this provides a completely independent test of the 
method of which did rely mainly on low energy information. Including more terms 
in the expansion, one obtains an improved approximation even in the low energy region. 
However, the quality of the "convergence" is significantly better for Ri and Rna than for 
Ra- Two reasons may be responsible for this difference: (i) In a high energy expansion it 
is presumably more difficult to approximate the 1/v Coulomb singularity in Ra than the 
mild Inf singularity in Rna and Ri. (ii) The function Ri can be approximated in the whole 
energy region 2m < a/s < 00 by an increasing number of terms with arbitrary accuracy. 
This is evident from the known analytical form of Ri, a consequence of the absence of 
thresholds above 2m in this piece. In contrast the functions Ra and Rna exhibit a four 
particle threshold at ^/s = 4m. The high energy expansion is, therefore, not expected 
to converge to the correct answer in the interval between 2m and 4m. For Rna this 
feature can be studied in more detail by separating Rna into the gluonic double-bubble 



terms in the C = 4 gauge |]18[ and a remainder. This separation is possible both for the 
semi-analytical result and the expansion (Fig. ^). The remainder approaches a constant 
both for a; — * and x — * 1. For < a; < 0.5 the agreement is perfect. It extends even up 
to a; ~ 0.9, a fact which is quite remarkable and surprising. 

For Rp the expansion is also shown in Fig. ^ Again one observes quick convergence 
for X between and 0.5. No (semi-)analytical result is available for the comparison 
with Rp. However, in the region below the four particle threshold an analytical result is 
available, based on the calculation of the form factor in 0. The four fermion contribution 
is expected to be small for energies just above 4m, corresponding to x just below x = 0.5. 



11 



0.1 0.2 0.3 0.4 0,5 0.6 0.7 0.8 0.9 1 0.1 0,2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 




Figure 2: The abelian contribution Ry , the non-abehan piece -R^^, the contri- 
bution from hght internal quark loops and the contribution Rp^ from the 
double-bubble diagram with the heavy fcrmion in both the inner and outer 
loop as functions of a; = 2m/ y^. Wide dots: no mass terms; dashed lines: 
including mass terms (m^/s)" up to n = 5; solid line: including mass terms 
up to (m^/s)^; narrow dots: semi-analytical result (except for Rp^). The scale 
//^ — has been adopted. 
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Figure 3: R\ ,i = A,NA,l,F for jj,'^ = s including successively higher orders 
in w? / s. The same notation as in Fig. || is adopted. The semi-analytical result 
is not shown. 
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Figure 4: — Rf\^ = 4) over x = 2m/y/s. Dotted: semi-analytical result; 
dashed: mass terms up to (m^/s)^; solid: mass terms up to (m^/s)^. 
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Figure 5: Comparison between the analytical result without 4-particle contri- 
bution {Rp{^^j.^, dotted line) and approximate result with terms up to order 
{m^/sY (dashed) and {w? / s)^ (solid line). 



Reasonable agreement between the two approaches is therefore expected in the region 
around x = 0.5. This is indeed observed in Fig. ^. 

5 Summary 

An algorithm has been described which produces the asymptotic expansion of the three- 
loop vacuum polarization automatically. It generates the relevant subdiagrams and as- 
signs the diagrams automatically to the programs MINCER and MATAD which evaluate 
the resulting massless propagator and massive tadpole integrals. The output is compared 
to the quartic terms obtained in and to the semi-analytical results of [||, confirming 
both these earlier results and the validity of the expansion down to fairly low energy 
values. 
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